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FOREWORD 


This  Final  Report  was  prepared  by  RCA  Laboratories,  Princeton,  New 
Jersey  under  Contract  No.  F30602-71-C-0356.  It  describes  work  performed  from 
9  June  1971  to  8  June  1972  in  the  Communications  Research  Laboratory, 

Dr.  K.  H.  Powers,  Director.  The  principal  investigator  and  project  engineer 
is  Dr.  D.  A.  deWolf. 

The  report  was  submitted  by  the  author  on  7  July  1972.  Submission  of 
this  report  does  not  constitute  Air  Force  approval  of  the  report's  findings  or 
conclusions.  It  is  submitted  only  for  the  exchange  and  stimulation  of  ideas. 

The  Air  Force  Program  Monitor  is  Lt.  Darryl  P.  Greenwood. 

PUBLICATION  REVIEW 

This  technical  report  has  been  reviewed  and  is  approved. 


SUMMARY  Am  ABSTRACT 


This  report  consists  of  two  parts.  The  first  is  a  review  and  discussion 
of  the  contents  of  the  three  Quarterly  Reports  under  the  contract 
effort  terminating  on  8  June  1972.  This  discussion  updates  that  work  by 
comparing  it  with  recent  other  work  and  new  data.  The  basic  results  deal 
with  the  fluctuating  properties  of  the  focal  area  of  a  focused  laser  beam 
in  turbulent  air,  viz., 

1.  The  average  focal  area  for  horizontal  propagation. 

2.  An  extension  of  (1)  to  slant -path  geometry. 

3.  The  power  spectrum  of  the  focal-spot  radius. 

The  second  part,  self-contained,  is  a  new  theory  of  plane-wave  irradiance 

scintillation.  It  contains  the  prediction  that  the  log-amplitude  variance 

,  -1/2  -1/6 

is  proportional  to  L  (C£)  in  the  saturation  regime.  The  theory 

is  an  extension  of  the  work  presented  recently  in  a  Special  Report  under 
the  present  contract  effort. 
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PART  I 


REVIEW  AND  DISCUSSION  OF  RESULTS 

The  results  obtained  and  reported  in  [l]  through  [4]  will  be  reviewed 
briefly,  and  discussed  in  context  with  the  most  recent  developments  else¬ 
where. 

A.  BROADENING  OF  A  FOCUSED  BEAM  IN  TURBULENT  AIR 

The  focal  spot  of  a  focused  laser  beam  scintillates  because  atmospheric 
turbulence  causes  irregularly  varying  undulations  in  the  rays  leaving  the 
optical  aperture  and  propagating  to  the  focus  (diffraction  limited  to  an 
area  characterized  by  a  typical  radial  length  ~  r_.  =  kL/ rQ  where  r^  is 
the  aperture  radius).  Three  specific  calculations  pertinent  to  this  effect 
have  been  given. 

1.  The  average  area  of  the  focal  spot  for  horizontal  propaga- 

2 

tion,  i.e.  uniform  C  . 

n 

2.  An  extension  of  (1)  for  slant-path  propagation,  using  the 

2 

altitude-dependent  model  of  Cn  (z)  summarized  by  Wyngaard, 
et  al[5] . 

3.  The  power  spectrum  of  the  focal-spot  radius. 

The  method  of  computing  (1)  requires  some  comment.  A  laser  beam  is  a  com¬ 
bination  of  modes  in  which  a  lowest-order  mode  usually  dominates [6] .  The 
ensuing  expression  for  its  electric  field  is  too  cumbersome  for  analytic 
work  and  consequently  an  approximate  approach  due  to  Schmeltzer[7]  is 


followed.  It  yields  a  result  for  the  irradiance[l]  that  contains  errors 
of  0(L/krQ2),  At  optical  frequencies  and  for  apertures  of  the  order  of 
1-meter  diameter  the  length  kr^  ~  10  km  and  therefore  this  approximation 
is  not  serious  for  applications  in  which  L  <  10  km  is  envisioned.  A  second 
choice  to  be  pondered  is  a  suitable  definition  of  focal-spot  area.  Several 
definitions  have  been  proposed: 

(a)  The  half-power  radius  r  '  defined  so  that  the  irradiance 

L 

at  (rL',L)  is.  related  to  that  on  the  main  axis  at  (0,L)  by 
<  I(rL'/L)  >  =  0.5  <  I(0,L)  >. 

II 

(b)  A  zero-th-moment  radius  r  , 

L 

(rLV  =  J*  d2p  <  I(p,L)  >/<  1(0, L)  > 

(c)  A  second-moment  radius,  r  , 

L 

\  H  J*  d2P  p2  <  KfU)  >/  J  d2p  <  i(p,l)  > 

(d)  An  irradiance-squared  definition,  r  (]T), 

ll 

rL2(I)  5  V  d2p  <  I<P^L)  >]2/  J  d2P  <  l2  (fU)  > 

Definitions  (a)  and  (b)  yield  rather  similar  results  which  require  a  num¬ 
ber  of  approximations  in  order  to  find  some  asymptotically  valid  analyti- 

y 

cal  expressions.  Definitions  (c)  and  (d)  can  be  modified  in  order  to 

2 

subtract  out  the  diffraction-limited  area  ~  r,  .  Definition  (d)  remains 

Lo  N 

intractable  analytically;  however,  Herrmann  and  Bradley[8J  have  used  this 


2 


definition  (without  the  averaging  procedure)  in  a  numerical  solution  of 

the  parabolic  wave  equation.  We  prefer  definition  (c)  because  it  yields 

[3] 

a  fundamental  result  that  can  be  interpreted  simply, 


r,  2  +  1.3  L3  K  ™  c  2. 
Lo  m  n  ' 


(1) 


where  K  is  the  microscale  wave  number.  This  result  has  been  compared  to 
those  obtained  from  (a)  in  [3],  and  it  was  concluded  that  the  difference 

t 

is  not  appreciable  in  spite  of  the  fact  (see  Eq.  (4)  of  [3])  that  r 

la 

2 

depends  upon  a  somewhat  different  combination  of  powers  of  ,  L,  and  k. 


Recently  [9]  data  on  beam  broadening  have  been  reported.  The  accompanying 

analysis  (and  presentation)  of  these  data  correspond  most  closely  with 

definition  (a).  Reasonable  agreement  between  optical-frequency  data  and 

calculations  exist,  but  the  infra-red  results  are  less  in  agreement.  As 

well  as  we  can  infer  at  this  time,  beam- spread  data  was  processed  around 

a  beam-wandering  center.  However,  the  analysis  based  on  expressions  (a) 

through  (c)  requires  that  no  distinction  be  made  between  beam  spread  and 

2  2  -1/3 

wander.  Specifically,  the  beam-wander  variance  <  tp  >  ~  Cq  Lb 

where  b  is  an  effective  beamwidth  radius.  The  measured  beam-spread  parame- 
2  2  2  2 

ter  is  <  6cp  >  =  <  9  >  -  <  cp  >  where  <0  >  is  the  angle  associated  with 
Equations  (3),  (4)  of  [ 3 ] .  Thus 

<  02>  ~  C  2  L  I  "1/3  for  0.4  k2  L  &  5/3  C  2  »  1 

no  on 

(3) 


C  12/5  16/5  k2/5 


for  0.4  k 


2 


L 


1 


3 


2 

We  think  the  trouble  may  lie  partly  in  the  fact  that  <  6cp  >  is  processed 
as  a  measure  of  beam  spread  and  then  compared  to  a  theoretical  expression 
for  <  02  >  as  in  (3). 

2  -13  -2/3 

First  consider  the  X  =  0.6-pm  case.  Let  &  =  1  cm.  C  «  10  m  ,  L  =  1  km. 

o  n 

2 

This  makes  Cn  L  correspond  to  the  measured  data.  With  these  values , 

2  5/3  2  2 

0.4  k  L  Z  C  ft*  1.  Hence  <  0  >  is  given  by  either  of  the  two  forms 
on 

-1/3  -1/3  -1/3 

in  (3);  they  are  numerically  the  same.  However ,  b  «*  LQ  s*  0.2  j6q 

2  2  -1/3  -1/3 

Thus  the  measured  quantity  <  (6cp)  >  ~  C  L  (j i  -  b  )~ 

9  -1/3  9  no  2 

0.8  Cr  L  ~  0.8  <0  >  is  interpreted  analytically  as  <  Q  >.  The 

error  is  at  worst  20  percent.  However,  for  X  =  10.6  pm,  we  find  that 

2  12/5  6/5  2/5 

<  0  >  is  now  given  by  Cn  L  k  ,  and  is  thus  reduced  at  most  by 

2/5  2  2 

the  factor  (0.6/10.6)  «  0.3.  In  that  case  <0  >  «  1.5  <  9  >  and  there 

2  2 

is  clearly  an  important  difference  between  <  6cp  >  and  <  0  >,  namely 
2  2 

<  6cp  >  «  0.3  <  0  >.  There  might  be  an  important  discrepancy  between 
measured  and  theoretical  beam  width  at  X  =  10.6  pm  for  this  reason. 

We  have  dealt  with  slant-path  propagation  in  [3].  The  basic  problem  is 

2 

that  the  knowledge  of  Cn  (z)  at  altitudes  z  above  several  hundred  meters 

ft  /n  /  /ft 

is  scant.  We  have  simply  extrapolated  the  z  and  z  models  of  [5]. 

9  2  2 

Instead  of  Eq.  (1),  we  obtain  for  rT  -  rT  =  rT  , 

L  LO  Lg 

2  2 

(s  / L)  Cn  (z)  downwards 

(2) 

(1  -  s /L)2  Cn2(z)  upwards 


9  9  1/3  L 

r_  z  «  3.9  L  K  f  ds 
Lb  m  J 


=  3.9  L2k  1/3  (*  ds 

m  J 


4 


Greenwood [10]  has  evaluated  these  integrals  in  terms  of  hyper geometric 

2 

functions.  In  [3]  we  present  some  asymptotic  limits  of  r  as  a  function 

of  the  altitude  Z  of  either  the  transmitter  or  receiver  (one  part  of  the 

2 

link  is  always  on  the  ground).  The  graphs  of  r  (Z)  versus  Z  must  be  used 
with  some  caution  when  Z  >  500  m,  but  they  serve  as  a  guideline. 


The  power-spectral  analysis  (3)  has  one  important  limitation;  it  is  re¬ 
stricted  to  the  angle  of  arrival  of  a  ray.  This  appears  to  us  to  be  cor¬ 
rect  for  a  focused  beam  as  explained  in  Section  IV. A  of  [3].  However,  some 
question  has  now  arisen  with  respect  to  the  reinterpretation  of  the  Gurvich, 
et  al.  data  presented  in  Figure  10  of  [11]  which  we  discussed  in  [3].  It 
is  possible  that  these  data  were  gathered  by  an  interferometric  technique, 
and  our  theory  is  not  applicable  to  an  interferometer  without  modifica¬ 
tion.  Nevertheless,  the  agreement  of  these  data  with  the  predictions  in 
[3]  appears  to  us  to  be  excellent,  so  that  possibly  our  reinterpretation 
remains  valid. 


B.  IRKADIANCE  SCINTILLATION  OF  A  PLANE  WAVE  IN  TURBULENT  AIR 

The  second  major  part  of  our  investigation  deals  with  what  is  perhaps  the 
most  fundamental  problem  in  the  field:  the  amplitude  fluctuations  o£  a 
plane  wave.  As  of  June  1971,  these  fluctuations  were  known  theoretically 
and  experimentally  to  be  governed  by  a  log-normal  probability  distribution 
determined  by  a  parameter  cr^  =  0.31  k?/6  ^1/6 

Cn^  in  the  homogeneously 

2 

turbulent  (horizontal  propagation)  case.  Under  conditions  in  which 


exceeds  unity  it  was  observed  that  the  log-amplitude  variance  saturates 

2 

then  decreases  as  increases.  No  satisfactory  theoretical  explanation 
has  been  given. 

We  have  summarized  extensively  in  a  special  report  [4]  what  seems  to  us 
the  most  promising  formalism  for  dealing  with  this  problem:  selective 
summation  of  terms  in  the  infinite  series  of  products  of  n-th  Born  terms 
in  each  moment  of  irradiance  I.  The  most  important  facet  of  this  report 
is  the  phys ical  interpretation  of  the  integrals  which  define  the  scatter¬ 
ing  processes;  i.e.  the  physical  interpretation  of  all  the  factors  in  each 
term.  Furthermore,  the  report  contains  error  bounds  that  appear  more 
accurate  than  found  elsewhere. 

This  work  has  led  to  an  approximate  solution  in  the  saturation  regime  of 
the  amplitude-fluctuation  problem:  the  amplitude  is  still  log-normal  and 
its  log-variance  varies  as  (c^2)-1/6.  The  major  break-through  was  to 
recognize  the  effect  of  ray-bending  upon  diffraction  effects:  a  prelimi¬ 
nary  version  of  the  new  theory  has  been  reported  in  [2],  The  full  theory-- 
which  includes  a  further  elaboration  of  the  work  in  [4]--is  appended  in 
this  report  in  technical-report  form. 

C.  ERRATA  TO  QUARTERLY  REPORT  NO.  3  (RADC-TR-72-119) [3] 

In  this  section  a  few  corrections  to,  and  additional  remarks  about  the 
results  presented  in  the  third  Quarterly  Report  [3]  will  be  added. 

(i)  The  numerical  coefficient  in  the  denominator  of  r  in  Eq.  (10) 

+1/2 

should  be  (1.3)  =1.14  (instead  of  1.4). 


6 


(id)  Figures  2  and  3  give  values  of  rQc  that  are  too  high  by  a  factor 
10.  To  correct  these  and  to  make  the  results  consistent  with 
Eq.  (10)  and  Figure  1, it  is  suggested  that  the  values  of  the  or¬ 
dinates  in  Figures  2  and  3  be  lowered  by  a  factor  10. 

(iii)  The  results  of  section  II  hold  for  air-to-ground  propagation,  where¬ 
as  those  of  section  III  hold  for  ground-to-air  propagation.  Thus, 
Figure  5  is  for  slant-path,  down  situations  and  Figure  6  (the 
caption  of  which  should  refer  to  laser  targets)  for  slant-path,  up 
situations.  All  text  in  sections  II  and  III  should  be  modified 

in  this  sense.  The  origin  of  the  original  confusion  is  the  re¬ 
ciprocity  invoked  in  interpreting  <  B(?,r^)B*(r ,r^) >  in  Eq. 

(A-3)  .  The  coordinate  s  in  Eqs.  (11)  and  (18)  should  have  its 

zero  at  the  receiver  end  because  of  reciprocity.  The  "dirty 

' 

windshield"  analogy  is  therefore  not  relevant. 

(iv)  An  «  sign  is  needed  in  Eq.  (36c). 

(v)  Table  II  is  for  UT  =  0  (not  UT»0) . 

(vi)  Instead  of  Eq.  (41b),  we  can  obtain  another  approximation  to  W0(a>) 
in  the  of2^3  region  from  Eq .  (34)  by  setting  q  =  5/6  and  1+x  «  x. 

It  yieldc 

n1/2  r(i/3)  2  2  _2/3 

Wa(co)  ~ - T7^ -  M(  1/ 6 , 1 ,  -&o  Mp  )*co 

9  (M)  ' 

where  M(=1F1)  is  a  Rummer  function.  This  formula  reduces  to  (36b) 
and  (38b)  for  oj  >  wT;  oj  >  Aw. 


7 


REFERENCES 


[1]  D.  A.  de  Wolf,  "Effects  of  Turbulence  Instabilities  on  Laser  Propaga¬ 
tion,"  RCA  Laboratories,  Quarterly  Report  to  Rome  Air  Development 
Center,  Griffiss  AFB,  New  York,  RADC-TR- 7 1-249  (October  1971). 

[2]  D.  A.de  Wolf,  "Effects  of  Turbulence  Instabilities  on  Laser  Propaga¬ 
tion,"  RCA  Laboratories,  Quarterly  Report  to  Rome  Air  Development 
Center,  Griffiss  AFB,  New  York,  RADC-TR-72-32  (January  1972). 

[3]  D.  A.  de  Wolf,  "Effects  of  Turbulence  Instabilities  on  Laser  Propaga¬ 
tion,"  RCA  Laboratories,  Quarterly  Report  to  Rome  Air  Development 
Center,  Griffiss  AFB,  New  York,  RADC-TR- 72-119  (April  1972). 

[4]  D.  A.de  Wolf,  "Strong  Amplitude  Fluctuations  of  Wave  Fields  Propagat¬ 
ing  Through  Turbulent  Media,"  RCA  Laboratories,  Technical  Report  to 
Rome  Air  Development  Center,  Griffiss  AFB,  New  York,  RADC-TR- 72-51 
(February  1972). 

[5]  J.  C.  Wyngaard,  Y.  Izumi,  and  S.  A.  Collins,  Jr.,  J.  Opt.  Soc.  Am. 

61,  1646  (1971). 

[6]  H.  Kogelnik  and  T.  I.  Lee,  Appl.  Optics  jj,  1550  (1966). 

[7]  R.  A.  Schmeltzer,  Quart.  Appl.  Math.  24,  339  (1967). 

[8]  J.  Herrmann  and  L.  C.  Bradley,  "Numerical  Calculation  of  Light  Propa¬ 
gation,"  Lincoln  Laboratory,  Laser  Technology  Program,  Report  LTP-10 
(12  July  1971). 

[9]  J.  A.  Dowling,  R.  W.  Harris,  M.  R.  Kruer,  and  T.  V.  Blane,  "Measure¬ 
ment  of  Wavelength  Effects  in  the  Propagation  of  Focused  Laser  Beams 
Through  Atmospheric  Turbulence,"  Paper  WEB  at  1972  Spring  Meeting  of 
O.S.A.  (11-13  April  1972). 

[10]  D.  P.  Greenwood, "Slant  Path  Spherical-Wave  Propagation,"  Paper  WE16 
at  1972  Spring  Meeting  of  O.S.A.  (11-13  April  1972). 

[11]  R.  S.  Lawrence  and  J.  W.  Strohbehn,  Proc.  IEEE  ,58,  1523  (1970). 


8 


PART  II 


STRONG  IRRADIANCE  FLUCTUATIONS  IN  TURBULENT  AIR:  PLANE  WAVES* 

By 

D.  A.  de  Wolf 
RCA  Laboratories 
Princeton,  New  Jersey  08540 
ABSTRACT 

Analytical  results  are  found  for  irradiance  statistics  of  a  plane 

wave  propagating  through  uniformly  turbulent  air.  They  arise  through 

a  new  error  analysis  of  approximations  made  in  selective  summation  of 

diagram  contributions  to  the  iterated  moment  equations  (in  integral  form) . 

The  results  are  determined  by  three  regimes  of  parameter  k^^L^^C  2 

n 

determined  by  its  location  with  respect  to  powers  of  the  micro-  and  macro- 

2  *1/6  2  a 

scale  Fresnel  numbers  ( Km  L/k)  ~  0,3  and  (kLQ  /L)  ~  104  (the  numer¬ 
ical  values  are  for  0.6-pm  radiation  and  a  horizontal  pathlength  L  ~  1  km) . 
The  modified-Rytov  result  is  found  in  the  lowest  regime.  Irradiance  I 
is  log-normal  in  the  intermediate  (saturation)  regime  and  the  loj -amplitude 

variance  decays  asymptotically  as  ( K  7^L3C  2)-1^6.  When  k7^6L*^6C  2 

m  n  n 

2 

far  exceeds  the  macroscale  Fresnel  number  kLQ  /L,  the  irradiance  tends  to 
an  exponential  distribution  (in  agreement  with  previous  results). 


* 

This  research  was  supported  by  the  Advanced  Research  Projects  Agency 
of  the  Department  of  Defense  and  was  monitored  by  Rome  Air  Development 
Center  under  Contract  No.  F30602-71-C-0356. 
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This  work  deals  with  the  statistics  of  the  irradiance  of  an  optical 
wave  propagating  through  uniformly  turbulent  air.  The  fundamental  sta¬ 
tistics  of  even  a  plane  wave  in  this  medium  have  been  predicted  only  for 
small  and  for  extremely  large  values  of  the  parameter 
ag  «  0.31  k  ^  L  ^  Cn  .  This  parameter,  a  function  of  pathlength  L, 

O 

wavenumber  k,  and  refractive-index  structure  constant  C  ,  appears  to  be  the 
one  that  governs  the  irradiance  statistics  for  a  plane  wave.  We  will 
present  new  results  for  the  hitherto  unsolved  intermediate  region  where 
ae  exceeds  unity  but  not  necessarily  by  very  much:  the  so-called  satura¬ 
tion  regime.  Also,  new  error  limits  will  be  given  that  enable  one  to 

separate  three  regimes  of  statistical  results  for  irradiance  I,  depending 

2 

upon  the  magnitude  of  and  the  micro-  and  macroscale  Fresnel  numbers. 

Tatarski  [l]  has  pointed  out  that  our  previous  result  [2]  yielding 

a  constant -plus-Rayleigh  distribution  for  I  is  incomplete  due  to  omission 

of  a  class  of  diagram  contributions.  We  have  found  this  criticism  to  be 

correct,  and  consequently  that  result  must  be  limited  to  the  radiowave 

case  [3]  of  large  macroscale  tresnel  numbers  L/kL^2 ,  in  which  case  the 

partial  summation  suffices  (but  a *  must  then  be  replaced  by 
2  5/3  2 

QL  0.2  k  LLq  Cn  ).  Subsequently,  we  amended  the  optical  case  [4] 
to  find  an  instantaneous  ray  solution  with  some  of  the  experimentally 
observed  characteristics  of  irradiance  statistics.  However,  the  solu¬ 
tion  only  predicts  the  statistics  of  I  globally,  and  the  ray  interpreta- 

2 

tLon  is  questionable  when  a  >1. 

e 

The  following  pertinent  developments  have  been  reported  in  this 
journal:  Sancer  and  Varvetsis  [5]  have  presented  an  irradiance  variance 
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vs  curve  that  saturates,  by  solving  an  approximated  equation.  Brown 
[6]  has  voiced  serious  objections  to  their  assumptions  and  approximations, 
and  their  results  appear  to  be  unconfirmed.  Young  [7]  attributes  satura¬ 
tion  to  beam  spreading  by  atmospheric  turbulence,  but  he  does  not  attempt 
to  give  a  full  analysis.  Molyneux  [8]  derives  equations  for  the  N-th 

it  2 

moment  of  I,  i.e.,  for  <i  >.  He  obtains  the  same  result  for  <  I  >  as 
in  [3]  and  therefore  Tatarski's  criticism  may  hold  here,  too.  A  related 
development  by  Brown  [9]  also  yields  general  moment  equations.  None  of  the 
above  works  allow  predictions  of  the  probability  distribution  of  I.  Furutsu 

[10]  derives  the  moment  equations  once  again  by  a  method  using  characteris¬ 

tic  functionals  of  the  refractive  index,  and  he  finds  that  the  logarithm  of 
the  irradiance  follows  a  Rice  distribution.  This  result  has  two  peculiari¬ 
ties:  (i)  the  irradiance  variance  goes  to  zero  as  the  gaussian-beam  field 

that  is  used  as  the  free-space  reference  field  approaches  a  plane  wave,  and 

(11)  the  exponent  of  the  mutual-coherence  function  is  assumed  to  be  a 
quadratic  function  of  the  separation  distance  of  the  correlated  observation 
points.  This  prediction  therefore  does  not  appear  to  apply  to  a  plane  wave 
in  turbulent  air. 

The  situation  in  the  USSR  is  somewhat  more  encouraging.  Important 
work  by  Tatarski,  Shishov,  Klyatskin,  and  others  prior  to  1971,  including 
numerical  solutions  [12]  of  the  moment  equations,  is  reviewed  by 
Barabanenkov ,  et  al. [11].  Since  then,  Gochelashvily  and  Shishov  [13] 

have  produced  further  analytical  results:  approximate  solution  of  the 

2  2 
<  I  >  equations  yields  the  Rytov  approximation  for  small  0g  ,  and 

-  asymptotically  -  an  exponential  distribution  of  I  at  large  distances. 
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Some  preparation  is  needed  to  present  the  result  of  this  paper. 
Consider  error  parameters  : 


■i2 ,  «mw'v 


“2 2  a  <L/H-o)ol2 


where  is  the  microscale  wavenumber.  For  kilometer  paths  at  optical 

-  .  2,  2  ? 
rrequencies,  o^  «  3^^  ,  and  the  numerical  factor  preceding  <j  varies 

slowly  with  L  and  k.  The  parametric  dependence  of  a  2  upon  L,  k,  and 

2  2  9  6 
cn  is  obtained  from  the  fact  that  a g  =  <  ^  >  where  X  is  the  real 

part  of  the  complex  phase  of  the  Rytov  approximation,  which  in  turn 
is  given  by 

=  t  dzjd£  (K,z)exp[-iK2(L-z)/2k]  ,  (2) 

Brt 

-2 

where  (2n)  dm  (K>z  )  is  a  Fourier-Stieltjes  increment  of  the  permittivity 
deviation  <5£  (p,z)  describing  the  Fourier  analysis  of  this  random  process 
in  any  z  =  constant  plane. 


Our  result  is: 

<  lN>  =  I0Nexp{2N(N-l)a£2[l-K)(a12)5/6]}  for  a*  «  1 

=  IoNexp[2N(N-l)ae2«0(ai2)-7/6]  for  «  1  «  (3) 

=  N'.  IoN[l  +  0 (a21/3e’a2)]  for  1  « 

The  electric  field  E  =  E^exp^j-a^2)  for  small  values  of  o^2.  Thus  the 
logarithm  of  the  irradiance  is  normally  distributed.  For  large  values 
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2 

of.  o 2  >  the  electric  field  is  governed  by  a  Rayleigh  distribution  of 

2  2 

amplitude.  There  is  an  intermediate  region  a ^  «  1  «  that  yields 

2 

a  log-normal  irradiance  asymptotically  as  0  such  that  the  log-amplitude 

2  -l/6 

variance  is  proportional  to  (ag  )  .  The  "saturation"  region  form  (for 

2 

1  £  Cj  )  lies  between  the  first  two  forms  of  (3)  and  it  appears  that  I 
is  log-normal  and  that  the  log-amplitude  variance  makes  a  transition 

from  a  linear  dependence  on  a  ^  to  a  dependence  upon  (a  in  agree - 

£  £ 

ment  with  experimental  evidence. 

The  result  is  obtained  by  the  method  of  selective  summation  of  dia¬ 
gram  contributions.  The  entire  Born  series  for  E  is  substituted  into 
N 

<  I  >  and  the  dominant  contributions  arising  after  a  diagram  renormal¬ 
ization  of  terms  are  summed.  The  error  analysis  is  therefore  a  crucial 
aspect  of  the  method.  A  geometrical -optics  interpretation  of  the 
change -of -angle  upon  scattering  allows  one  to  estimate  the  error  and 
to  make  very  significant  approximations. 


I.  THE  BORN  SERIES  IN  SMALL-ANGLE  SCATTERING  APPROXIMATION 

A  very  brief  derivation  of  the  parabolic  wave  equation  in  the  form 
most  useful  for  selective-summation  methods  will  be  given  because  the 
source  of  the  approximations  needs  to  be  specified  carefully.  The  basic 
equation  for  a  plane  wave  propagating  from  the  origin  plane  z  =  0  to 
the  point  of  observation  r  =  (0,0, L)  is  [4]: 

2 

B(r)  =  1  +  ~  /d3rltfd6  (K1,z1)G(r-r1)exp(-iK1<P1)B(r1)  (4) 
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»here  coordinate  ^  d£  (Kj.z)  is  the  two-dimensional  Fourier- 

Stieltjes  Increment  when  the  random  variable  6£  (i)  is  Fourier-analyzed 
in  a  z  -  constant  plane,  B(r)  is  the  electric  Field  E(r)  divided  by  the 
equivalent  free-space  Field  E  <r)  =  eXp(ikz),  and  G(r-r.)  is  the  Green’s 
function  multiplied  by  the  factor  <exp[-ik(Z-2l)]} .  The  d2^  integration 
in  (4)  can  be  performed.  Let  us  write  B^)  =  AOr^expficp^)],  „ith 
real  functions  A  and  <p.  Then  stationary-phase  analysis  [4]  shows  that 
the  only  region  in  the  ^  =  const .  plane  of  importance  is  the  immediate 
vicinity  of  the  stationary  points ^  that  obey  the  equation 


(p-pp/lr-ril  -  -K./k  +  k'V^rj) 


(5) 


where  VT  is  the  transverse  derivative,  in  this  case  at  r 


However,  the  d^  integration  is  not  carried  out  for  (4),  but  for 
_ _ _  00 

Its  iterated  form  B  -  Z  B^  and  it  cau  be  done  without  stationary- 
phase  analysis  to  obtain, 

B”  ‘  A  d2J  fa  <Vj> 

X  Jl<l'Q">2/k2>1/2e,!p(lkf(1'Qm2/k2)1/2|a2ml-atm])  .  (61 


X 


where  Azm  ZQ  =  L,  andjgm  is  a  shorthand  notation  that  does 

not  show  the  dependence  on  index  n.  If,  however,  a  stationary-phase 
analysis  is  done  on  each  d2pm  factor  of  Bn  obtained  from  (4)  by  the  n-th 
iteration,  even  though  we  did  not  use  this  analysis  in  deriving  (6), 
we  find  that  the  only  regions  in  the  z n  =  const,  planes  of  importance 
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(7) 


are  the  immediate  vicinities  of  the  stationary  points  p  that  obey, 

~m 


(o  i“Pm)/|r 


i  "r 

m-1  -m 1 


There  is  an  important  connection  between  (7)  and  (5)  that  we  wish  to  ex¬ 
ploit  in  this  work.  However,  first  we  discuss  two  approximations  always 
made  in  (6): 


Approximation  (i):  We  set  z  <  z  ,  for  all  z  ,  i.e.,  backscatter  con- 

m  m-1  m 

tributions  are  ignored.  The  effect  has  been  discussed  before  by  many 

workers,  e.g.  in  [ 1 1 ] .  It  can  be  understood  by  regarding  the  spectrum 

4>(K,k  )  of  6£  (p,z).  The  ratio  of  the  variances  of  the  d£(K.,z.)  fac- 
—  z  —  —J  J 

tor  in  (6)  when  dz^  is  integrated  over  a  distance  -Az  and  H-Az  is  propor¬ 
tional  to  the  ratio  $(K,2k)/®(K,0)  where  2nL  *  <  K  <  K  .  In  optics 

•a*  o  m 

flCK^k)  «  0  for  all  practical  purposes  and  the  dz^  integration  over  -Az 
is  negligible. 


2  2 

Approximation  (ii):  Terms  of  0(Qm  /k  )  are  ignored  compared  to  unity. 

2  2  2  2 

When  n-m  «  10  one  might  estimate  Q  /k  £  k  /k  «  1.  Unfortunately, 

m  m 

2  2 

large-n  terms  are  sometimes  required,  n-m  can  become  large,  and  Q  /k 

2  2 

may  greatly  exceed  K 1  /k  .  Furthermore,  there  is  a  cumulative  effect  of 

22  22  1/2 

errors  due  to  dropping  all  Qm  /k  in  (1-Qm  /k  )  1  in  (6). 


2  2 

A  great  aid  in  understanding  the  effect  of  terms  of  0(Qm  /k  )  com¬ 
pared  to  unity  in  (6)  is  the  geometrical  interpretation  of  the  vectors 

Q  /k  given  by  (7).  It  states  that  |Q  /k I  is  the  (sine  of  the)  angle  of 
«*m  ">m 

the  stationary-phase  ray  determined  by  the  increments  dg  (K  ,z  )***d£(K  ,z  ) 

—m  m  n  n 

for  given  zm>  Let  us  consider  (7)  for  m=l,  and  simultaneously  for  all 

n  >  1.  Alsc,  all  dg  (K. ,z.)  are  allowed  to  vary  in  all  terms  for  j  >  1 
~  -'J  J 
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(i.e.,  the  j  =  1  increment  is  not  varied).  The  total  effect  of  (p-p, )/ |r-r. [ 

of  (7)  in  all  these  terms  must  be  identical  to  that  in  (5).  The  total 

collection  of  allj^/k  (with  correct  sign  preceding  it),  which  we  denote 

by  (Q./k)  must  then  be  identical  in  effect,  as  all  d&  (K.,)z.)  are  are  varied, 

.  •**'  J  J 

to  the  right-hand  side  of  (5).  Because  the  K./k  term  subtracts  out,  we 
have, 

(22/k)  =  k’^cpCrj)  (8) 

If  we  iterate  (4)  m-1  times,  then  find  the  stationary -phase  equivalent 
of  (5)  for  in  the  last  term,  and  then  apply  the  above  comparison 

with  (7)  in  similar  fashion,  we  obtain 

LW11  '  k'Vrn,>  <9> 

in  the  understanding  that  d£  (K^,z^)  is  varied  for  all  j  >  m,  and  that  all 
terms  Bn  with  n  >  m  are  considered  simultaneously.  This  is  an  extremely 
valuable  identification  because  the  right-hand  side  of  (9)  can  be  iden¬ 
tified  with  the  angular  deviation  of  a  ray  in  a  geometrical -optics 
medium: 

k‘V(rm)  =  d»(zm)/dzm  <10> 

where  the  (for  small  angles  quite  permissible)  approximation 
dp/dr  »  dp/dz  has  been  included.  Note  that  we  cannot  replace  any  par¬ 
ticular  Q  . . /k  by  dp/dz  in  (6).  But  we  can  estimate  the  effect  of  the 
ntTL 

(Q./k)  in  (6)  by  regarding  the  behavior  of  dp(z  )/dz  through  applica- 
vnm  m  m 

tion  of  (9)  and  (10).  No  matter  how  complicated  9(rm)  actually  is,  the 
phase  derivative  is  determined  by  the  directions  of  the  rays  through 
the  point  r  .  Thus 
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(Ha) 


2 

dp(z  )/<iz  =  “  V_,  /  m  dz  6e  [p(z),z] 


after  making  the  small-angle  approximation  of  straightening  rays.  The 
right-hand  side  of  (11a)  is  a  statistical  quantity  that  approaches  a 
normal  variable  with  zero  mean  and  with  variance, 


<  (dp  /dz  )  > 
— m  m 


K  1/3z  C  2 
m  m  n 


(lib) 


in  a  turbulent  medium  with  the  usual  modified  Kolmogorov  spectrum 

<  (5  e)2>«(K,0)  ■  C  2K-11/3exp(-K2/K  2).  When  K  1/3LC  2  «  1,  we  can 

n  m  m  n 

2  2 

be  sure  that  all  Q  / k  are  small  to  this  order  because  the  statistical, 
m  ’ 

nearly  normal,  variable  |dp/dz|  will  seldomly  exceed  several  times  its 
ms  value  determined  by  (lib).  Thus,  the  approximation  amounts  to  the 
neglect  of  diffusion  effects:  the  cumulative  scattering  at  large  angles 
is  ignored.  In  short,  approximation  (ii)  is  based  on  the  statistical 
(geometrical -optics)  estimate 


<,«>2  s  ocC'H2)  • 


m 


(12) 


1/3  2 

and  when  K  L  C  «  1  we  find  that  (6)  reduces  to, 
m  n 

Bn  ‘  jll  ^2  /  3  dzj/dS  -S=l  exI,[-1<!n12(zm-rzn,)/2kl 


X  [l+0(/C  1/3l  c  2)] 

m  n 


(13) 


There  is  an  extra  condition  hidden  in  (13).  When  the  factors 

2  2  1/2 

(l"Qm  /k  )  are  expanded  in  a  Taylor  series  in  the  exponential  in  (6), 
the  leading  terms  are  estimated  from  ikzm[(l-Qm2/k2)1^2-(l-Qm+12/k2)1^2] . 


m 


2  2  3 

The  first  one  is  given  in  (13),  and  the  second  one  is  0(K  Q  ,,  L/k  ) 

m  m+1 
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and  -  in  accordance  with  all  sagittal  approximations  -  it  must  be  much 
less  than  it .  Using  (12),  we  find  this  error  to  be 

°(K  2Q.  1 2Ii/k3 )  <  0(L2k  1/3k_1L  C  2) 
m  m+1  m  on 

2  1/3  -12 

hence  L/C  k  L  C  «  1  for  validity  of  (13).  However,  because 
m  o  n  J 

2  1/3  2 

L  «  kL  we  note  that  this  condition  is  implied  by  K  LC  «  1, 
o  r  m  n 

and  therefore  it  need  not  be  stressed. 

Now  we  return  to  (13),  which  is  equivalent  to  the  parabolic  equa¬ 
tion  [11],  and  rewrite  it  as 

Bn  *  d'j/ d£ 

n  (Wa) 

X  L5-,  exp  [  -iK  2  (L -z  ) / 2k ] exp [ -iK  *Q  . .  (L -z  ) /k ] } 
in— i  m  m  ^  m  i^ntTi  m 

The  two  exponential  factors  can  be  ignored  in  the  high-frequency  limit 
k  -»  “.  As  we  have  shown  previously  [3],  and  as  Modesitt  [14]  has  recently 
pointed  out  again,  (14)  can  be  summed  to  yield  the  WKB-phase  integral 
in  the  form  also  known  as  the  Moliere  approximation.  The  Rytov  approxi¬ 
mation  follows  if  only  the  second  exponential  factor  for  each  m  is  ignored. 

2 

The  modified -Rytov  approximation  given  in  (3)  for  «  1  follows  by 

N 

essentially  ignoring  this  exponential  factor  in  <  I  >  for  N  >  1  rather 
than  in  B.  In  fact,  the  study  that  follows  is  centered  around  the  role 
of  these  exponential  factors  when  (14)  is  substituted  into  <  IN>  . 

Let  us  write  (14a)  formally  in  a  short -hand  notation: 
n  z 

B  =  II  f m"1  dz  d£  (m)F  (m)  (14b) 

n  m=o  J  m  n  v  ' 

o 
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I 


Let  us  first  consider  <  IN>  for  N  =  1.  It  consists  of  a  sum  of  terms 


<  B  B  *  >  in  which  p  and  q  are  summed  from  zero  to  infinity.  Let  us 

p  q 


consider  a  particular  term  for  which  p  +  q  =  n.  We  write  this  term,  using 


(14b),  as 


f  dZl...  Jp_1  dzp  <V  •  •  /,'1  «  (l)Fp(l)...d£(q)Fq*(q) 


For  given  z^”’z  .z^  •  • -z  1  we  can  rearrange  the  factors  to  become 


/v-/ 


dz  da  (1")F  (1")  ...da  (n")F  (n") 


where  z  "•■•z  "  is  a  permutation  of  z  •••z  ,z.’'**z  1  that  leaves  the 

In  P  1  q 

ordering  z.  ,  >  z.  and  z.  .’  >  z\  unchanged.  Likewise  P  (i")  is  one  of 
°  1-1  L  J  L1 


the  F  or  F  *  factors  depending  upon  which  coordinate  is  permuted  into 

p  q 


z.".  There  are  exactly  n'./p'.q!  integrals  of  type  (16)  that  constitute 
(15).  A  similar  procedure  can  be  carried  out  for  N  >  1  but  obviously 
it  is  much  more  cumbersome  to  describe.  The  point  is  that  (16)  is  now 


an  ordered  integral 


Two  simplifications  of  the  n-point  integral  in  all  terms  (16)  are 


possible  regardless  of  the  joint  probability  density  of  da  (l)...da(n) 


provided  the  binary  correlation  length  2  «  L  (the  length  2  is  of  the 


order  of  the  macroscale  L  )  . 

o 


(i)  The  n-point  correlation  can  be  separated  into  n'./2n  ^(n/2)'.  products 


of  binary  correlations  <  da  (i)da(j)  >  .  All  products  containing 


2  2 

higher -order  correlations  are  negligible  by  errors  0(1^  /L  )  to 


2  2  2  2 

some  power,  and  cumulatively  by  0(k  Lq  e  )  where  e  is  a  short -hand 


- - 


notation  for  the  permittivity  variance.  The  proof  is  given  in 
[3]  and  more  exhaustively  by  Brown  [15]. 

(ii)  Of  the  remaining  products  of  binary  correlations,  only 

<  d£  (l")d&(2")>  •••<d£  (n-l")  d£  (n")>  is  significant.  The 

others  contain  at  least  two  extra  restrictions  of  intervals  Az 

2  2 

to  a  distance  ~  L  which  introduces  reduction  factors  of  0(L  / L  ) 

o  o 

2  2  2 

per  product  and  0(k  Lq  e  )  cumulatively  again.  Furthermore  there 

2  2 

are  some  errors  0(Q^  /k  )  due  to  the  factors.  These  are  negligible 
by  virtue  of  (9)-(ll). 

In  summary,  the  normalized  electric  field  B(=E/Eo)  can  be  expressed 

as  a  sum  of  Br  over  all  integer  n,  which  are  given  by  (14).  These  are 

inserted  into  <  1^  >  =  I  (BB*)^>  and  the  above  properties  (i)  and 

(ii)  are  applied  to  the  many  ordered-integral  terms  which  go  into  this 

2  2  2 

N-th  moment.  The  results  are  then  restricted  by  errors  0(k  L  e  )  and 

o 

1/3  2 

°(^m  LCn  )•  Figure  1  is  very  useful  for  illustrating  these.  The  cross- 

hatched  areas  are  those  where  these  errors  are  large.  The  other  regions 
2  2 

of  this  e  (or  Cn  )  vs  L/Lo  plot  are  useful  for  defining  parameter  regimes. 

Optical  propagation  appears  to  be  restricted  to  regions  I  and  III.  The 

theoretical  difficulties  lie  in  the  transition  from  III  to  III  ,  because 

w  s 

2 

this  is  where  increases  above  unity.  The  other  areas  are  more  ger¬ 
mane  to  radiowave  propagation  and  they  need  not  be  discussed  at  present. 
Previous  solutions  for  <  IN  >  in  [2]  and  [3]  are  restricted  to  regions 
I,  II,  and  IV  (when  N  >  1) . 
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- (L/K£?)'/S  a ,2  =1 

K?/6  Ll,/6Cn2 


Figure  1.  Parameter-region  diagram  for  approximate  solutions  to  the  wave 
equation  in  turbulent  air.  Cross-hatched  regions  indicate 
breakdown  of  the  parabolic  equation  and/or  of  important  statis¬ 
tical  assumptions. 
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II.  THE  MODIFIED-RYTOV  APPROXIMATION 
2 

Let  us  consider  <  I  >  .  It  is  determined  by  the  sum  over  all  in¬ 
teger  n,  m,  p,  q  of  <  B^Bm*BpB^*>  •  Let  n+m+p+q  =  2M.  The  statistical 

simplification  explained  in  connection  with  (15)  allows  one  to  represent 
2 

<1  >  by  a  sum  of  diagram  contributions,  each  diagram  corresponding 

not  to  n,  m,  p,  and  q,  but  to  the  half  sum  M.  Each  diagram  consists  of 

horizontal  z-axes  placed  under  each  other,  one  for  each  B  factor  (dashed 

lines  for  B*  axes)  Each  two-point  correlation  <d£  (i)  d£  (i+l)>  is 

represented  by  a  "bead"  on  an  axis,  or  by  a  "rung"  between  two  axes, 

depending  upon  which  B-factor  the  random  increment  d£  (or  perhaps  its 

complex  conjugate)  is  part  of.  Figure  2  gives  some  examples.  It  can  be 
2 

seen  that  <  I  >  is  given  by  the  sum  of  all  topologically  different  dia¬ 
grams  with  M  features  (M  -  0,  1,  •••).  Note  that  an  M  =  2  diagram  with 
a  bead  to  the  left  of  a  rung  is  topologically  different  from  one  with  a 
bead  to  the  right  of  a  rung.  The  diagram  rules  have  teen  explained  be¬ 
fore  [2,3]. 


Consider  an  M-feature  diagram  and  keep  features  2,  3,  •••,  M  fixed. 

The  first,  feature  can  then  be  one  of  the  ten  shown  in  Figure  3.  We  now 

utilize  the  following  property  in  <  B  B  *B  B  *>  with  (14) 

n  m  p  q  ■ 

<dS  (1)  d£  «)>  -  fcr2e2* ,<J ,)d2K.d2K.  6,<K.  +  K.),  (17) 

The  minus  sign  in  the  (two-dimensional)  Dirac  delta  factor  holds  if  i 
belongs  to  a  B  and  j  to  a  B*  factor,  or  vice-versa,  and  the  plus  sign 
in  the  two  other  cases.  The  function  ^(K,^)  can  be  Fourier  transformed 
wfth  respect  to  Az  into  the  well-known  three-dimensional  spatial  turbulence 
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(a)  A  BnBm  RUNG  AT  Zp 


Zp+2 

(b)  A  B„Bm  RUNG  AT  Zp 


Figure  2.  Two  examples  of  diagrams  symbolizing 
contributions  to  <  I 2  > . 


Figure  3.  Ten  possible  locations  of  a  feature 
of  an  <  1 2  >  diagram. 


spectrum  $(K,kz).  In  the  present  work,  the  only  form  of  the  spectrum 
that  will  appear  is  $(K,G)  which  is  written  as  ®(K)  for  convenience. 
The  first  factor  of  the  sum  of  ten  M-feature  diagrams  thus  obtained  is 


(ke/4it)2  JL  dzJdh^Kj 

X  {[C1(12)+C1(34)+C1(14)+C1(23)-2]  -  [Cx ( 13 )-K31  (24) ]C  ) 
C  s  cos[K12(L-z)/k] 

C^ij)  =  cos[K1*(Q2(i):§2(j))(L-z1)/k] 


(18a) 


The  vectors^Q^^  1°  the  definitions  of  C^(ij)  pertain  to  the  sum  of 

JKp-vectors,  p  >  2,  connecting  the  i-tU  B  axis  of  the  diagram  to  all 

other  axes.  Note  in  particular  thatj)2^  andj^^  are  kept  fixed  for 

all  ten  diagrams  generated  by  allowing  the  m  =  1  feature  to  take  ail 

ten  choices  of  Fig.  3.  Note  also  that  these  ten  diagrams  have  already 

incorporated  some  symmetry  features  with  other  sets  of  ten  diagrams  to 

yield  only  cosine  factors  in  (18a).  For  example,  the  contribution 

<  BmBn*B,Ai*  >  can  be  taken  together  with  the<B  B  *B  B  *  >  contribution. 
* “  11  P  H  n  m  q  p 

The  second  and  other  factors  following  (18a)  cannot  be  written  out  in 


similar  fashion  because  (18a)  contains  vectors  KL,#,'K  .  If  we  were 

— -m 

to  change,  say,  the  second  feature  (m  =  2),  then  this  would  modify  some 
of  the  Cj(ij)  factors  in  (18a).  However,  we  can  also  take  a  particular 
M-diagram  without  allowing  any  of  the  features  to  be  permuted,  and  then 
we  obtain  a  contribution, 


/d2^  •  •  •  fd 


M 


.2  2 

k  e 


,|CM  mSi  2 

16jt 


2 

/m_1  dzm<t(Km)  X  (-l)i+j_1Fm(ij)  (18b) 
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where  F^(ij)  stands  for  one  of  the  cosine  factors  in  {•••}  of  (1.8a), 
but  for  any  m,  when  i  3*  j,  and  Fm(ii)  =  1/2.  All  the  bookkeeping  diffi¬ 
culties  lie  in  the  Cm< i j )  factors.  We  will  use  both  forms  (18)  in  the 

following.  Note  finally  that  (18a)  corresponds  to  the  moment  equation 
2 

for  <  I  >  preferred  by  other  workers. 


The  modified-Rytov  approximation  is  obtained  by  setting  all  C  (ij) 
Consider  m  =  1  first  and  set  C^(ij)  =  1  in  (18a),  which  then  reduces  to 


=  1 


.2  2 

k  e 


“  w  pq 

J  d*!  J  dK^K^d-cosCK^a-z^/k]) 


(19) 


This  reduced  form  of  (18a)  is  independent  of  vectors  K_  • • *K  .  Therefore 
the  m  =  2  feature  can  be  treated  in  exactly  the  same  way  as  the  m  =  1 
feature  in  (18a)  ;  after  permutation  of  the  second  feature  a  common  fac¬ 
tor  identical  to  (18a)  results  except  that  we  have  C2(ij)  factors.  Now 
we  set  C2(ij)  =  1  to  obtain  a  form  similar  to  (19).  The  product  of  all 
these  factors  for  m  =  1  to  m  =  M  must  exhaust  all  possible  M-feature 
diagrams .  One  need  only  sum  from  M  =  0  to  M  =  »  to  obtain  <  I2  > 

1  *6 « 


<i2> 


m  m  m 


m 


'  "I  Tr/"'1  /<K„M(KJ(1-C_> 

,  ,2  2 

_  t  2  rk  e 

=  I  ex  - 

o 


P[^T  £  dz  /  «1K1*(K1)(l-C1)]  (20) 


=  I  2exp(4a  2)  . 
o  e 


The  procedure  is  easily  generalized  for  <  IN  >  with  N  >  2.  An  expres¬ 
sion  similar  to  (18a)  can  be  written  down;  the  main  difference  is  that 
^(l2)^"  "b  C^(23)  is  replaced  by  a  sum  of  N  terms  corresponding  to 
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<  BB*  >  rungs,  the  numerical  term  -2  is  replaced  by  -N,  and  ^(13)4^(24) 
is  replaced  by  N(N-l)  terms  exhausting  all  <  BB>  and  <  B*B*  >  rungs. 

It  need  not  be  written  down.  The  same  procedure  is  easily  observed  to 
yield, 

<  lN  >  =  tcN  exp[2N(N-Dce2]  .  (21) 

In  this  fashion,  the  first  of  Eqs.  (3)  is  easily  derived.  The  remain¬ 
ing  analysis  of  this  result  is  for  the  purpose  of  estimating  the  error 
induced  by  setting  all  Cm(ij)  =1.  We  shall  do  this  systematically, 
first  for  m  =  1  then  for  m  =  2 ,3, • • • ,  etc.  It  is  preferable,  however, 
to  simplify  (18a)  first  by  means  of  an  intermediate  renormalization  which 
simplifies  the  diagram  bookkeeping. 

RENORMALIZED  BEADLESS  DIAGRAMS 

There  is  a  contribution  of  a  numerical  term  -2  in  the.  {  }  factor 

in  (18a)  due  to  the  four  possible  positions  of  beads  in  the  first  factor. 
It  is  therefore  quite  easy  to  renormalize  the  integrals  by  summing  out 
the  contribution  of  all  beads,  and  the  advantage  of  the  new  expressions 
is  that  the  renormalized  diagrams  have  rungs  only.  Let  us  start  with 
any  diagram  without  beads  interpreted  by  meant',  of  (18b).  Consider  in 
this  diagram  any  axis  segment  between  two  contiguous  rungs.  Now  generate 
a  new  set  of  diagrams  by  placing  zero,  one,  two,  —  beads  on  this  par¬ 
ticular  segment.  The  sum  of  the  contributions  from  these  new  diagrams 
can  be  seen  from  (18b)  to  correspond  to  an  extra  factor 

k2  2 
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in  the  contribution  (18b).  We  can  do  the  same  thing  for  every  axis  seg¬ 
ment  of  the  original  diagram.  The  result  is, 

(i)  we  have  generated  all  possible  diagrams  with  M  rungs  and  any 
number  of  beads; 

(ii)  we  have  generated  a  factor  exp(<i!L)  for  every  horizontal  axis, 

2 

i.e.,  a  factor  exp(-4oL)  for  every  <  I  >  renormalized  dia¬ 
grams.  (The  explicit  form  of  (XL  is  given  in  the  introduction.) 

Clearly  there  is  nothing  special  about  M,  and  the  net  result  is  that 
we  can  restrict  ourselves  to  beadless  diagrams.  An  M-rung  renormalized 
diagram  is  interpreted  as  follows.  The  corollary  of  (18a)  is 


exp(-4d!L)  *(k€/4n)^  dZjJd^K^(K^) 

X  {[C1(12)-K31(34)-K31(14)-K31(23)-C1(13)-C1(24)]  (2 

+[c1(13)-k:1(24)](i-c1)}  , 

and  the  corollary  of  (18b)  will  not  be  written  down  here  because  it 
differs  only  by  a  factor  exp(-4o!L).  Note  that  we  have  reordered  the 


terms  somewhat  in  (22),  aside  from  the  renormalization.  Steps  (19)-(21) 

are  easily  repeated,  using  (22).  The  only  difference  is  that  {•••}  re- 

2 

duces  to  2  +  2(1-Cj'  ■  -ilch  yields  a  factor  4o!L  +  4 0£  at  all  places  in 

2 

(20)  and  (21)  where  previously  only  4 o£  appeared.  This  cancels  the 
factor  exp(-4o!i).  There  is  a  definite  advantage  to  the  reordering  in 
(22)  which  will  become  apparent  in  the  error  analysis.  In  order  to  in¬ 


vestigate  the  effect  of  setting  C^(ij)  =  1  in  (22),  we  utilize  (23a)  as 
the  prototype  for  analyzing  any  of  the  terms  in  the  second  row  of  (22), 


and  (23b)  for  the  third  row  of  (22). 


(ke/4 nt  )2  JL  dzJ*dZK®(K)cos  [K-/5Q(L-z )k] 

(ke/4n)2  fL  dzJd2K®(K)(l-cos[K2(L-z)/k]  }cos[K;^(L-z)/k] 


(23a) 

(23b) 


o 

Here,  iQ  is  a  short-hand  notation  for  The  error  analysis 

W"  IMr  Z 

is  not  affected  if  (L-z)  in  (23)  is  replaced  by  L,  to  obtain 


(k2e2L/8rt)  f0  dKK4(K)JQ(K4QL/k) 


(ke2  L/8jt)  J°°  dKK® (K ) Jq  ( K^QL/k )  ( 1  - ( k/ K  L ) s  in ( K  L/k  ) )  , 


(24a) 

(24b) 


where  J  is  a  Bessel  function.  From  here  on,  the  Kolmogorov  spectrum 
o 

[16]  is  used,  with  the  von  Karman  modification  for  wavenumbers  K  <  Lq 


®(K)  «  15.7Lo3(l+K2Lo2)“U/6exp(-K2/Km2) 


(25) 


The  numerical  factor  in  (25)  differs  from  that  used  by  Tatarski  [16] 
due  to  a  different  Fourier -transform  normalization.  Insertion  of  (25) 
into  (24)  yields  forms  which  serve  as  the  basis  for  error  analysis: 

(15.7/4It)k2e2LLo[o25/6K5/6(a2)/211/6r(ll/6)]  ,  (26a) 

with  a0  =  £QL/kL  , 
z  o 

4o  2b  (a.)  with  a  =  £Q(L/k)1/2.  (26b) 

6X1  L 

Here,  K5/g^a2^  is  a  Bessel  function,  T(ll/6)  is  a  gamma 

function,  and  bx(a^)  is  related  to  a  confluent  hypergeometric  function, 
specifically  by  Eq.  (47.33)  of  [16].  Only  the  following  two  limiting 
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forms  of  this  function  are  required  in  this  work: 


b^)  =  1  -  2.36(c1V/6  as  ot  -  0 
=  -0.0242(o12/4)"7/6  as  Oj 


(27) 


There  is  a  small  detail  in  the  ->  0  form  regarding  the  magnitude  of 

9  2 

0^  with  respect  to  the  inverse  microscale  Fresnel  number  /L  that 

is  of  no  importance  here  (it  is  assumed  that  does  not  decrease  below 
this  number) .  With  (26)  it  is  easily  seen  for  small  £Q/k  that  (26)  re¬ 
duces  to, 


4c*L[l  +  0(a22)5/6] 

(28a) 

4ae2[l  +  0(a12)5/6] 

(28b) 

A  factor  2  has  been  multiplied  into  (26)  and  (28)  to  make  (24)  correspond 
to  the  sum  of  terms  in  (22).  Furthermore,  the  definitions  imply  that 
a ^  «  Oj2  so  that  it  suffices  to  require  a^2  «  1  in  (28b)  in  order 

that  the  modified-Rytov  expressions  (20)  and  (21)  follow. 

2 

Wich  the  development  of  (9) -(11) ,  it  then  follows  that  can  be 

identified  with  the  parameter  of  Eq .  (1).  There  is  one  difference.  In 

considering  (22)  for  (^(ij),  it  is  necessary  to  consider  simultaneously 

all  diagrams  with  M  >  1.  In  that  casej^/k  can  be  identified  again  with 

k_1Vmcp(r.),  and  <2Q„/k  is  then  certainly  less  than  2k  *V  cp(r-);  it  can 

even  be  zero  if  the  two  runged  axes  to  which  belong  form  an  uncon- 

2 

nected  subdiagram.  However,  in  the  worst  case,  we  can  replace  (<5Q/k) 

by  4Lk  ^3C  2  (ignoring  a  factor  of  order  unity)  and  thus  convert  the 

2  2 

above  defined  parameters  and  o2  to, 
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(29) 


„,2«.  kL2*  1/3C  2  =  (k  2L/k)1/6(k7/V1/6C  2) 

I  m  n  m  n 


022«-»(L/kLo2)012  . 


III.  THE  SATURATION  REGIME:  1  <  k7^6LU^6C  2  «  kL  2/L 

n  o 

The  results  of  the  previous  section  indicate  that  the  log-amplitude 

(or  log -ir radiance)  variance  is  a  linear  function  of  the  parameter 

a  2  ~  k^^L^^C  2.  All  measurements  to  date  show  a  variation  from  this 
e  n 

2 

behavior  when  a  >  1:  a  saturation  effect  occurs  in  the  sense  that  the 
e  * 

2 

variance  attains  a  maximum  value  close  to  o£  =1.  We  shall  define  a 

2  2 

saturation  regime  by  the  inequality  1  <  0£  «  kLQ  /L,  and  it  corresponds 

2  2 

reasonably  well  to  the  slightly  sharper  inequality  a2  «  1  «  a ^ 
whenever  the  identification  (29)  is  utilized. 

2  2 

Clearly  o  «  kL  /L  does  imply  that  £Q/k  in  the  first  set  of 

€  O 

2  S/6 

terms  of  { - )  in  (22)  is  small  enough  to  set  C^(ij)  =  1  +  0(a2  ) 

for  the  same  reasons  as  set  forth  in  Eqs.  (23) -(28).  Consequently, 

(22)  reduces  to, 

exp(-4o(L)'[(kt74n)2o/i  dz1/d2K1®(K1){2  +  [^(13)4^(24)] [l-C^})  (30) 

Through  an  obvious  symmetry  property  that  maps  every  diagram  of 
<  BB*BB*  >  into  one  of  <  B*BB*B  >  ,  and  vice-versa,  C^(24)  may  be  replaced 
by  C^(13).  The  filter  factor  {— }  in  (39)  reduces  to 
2{1  +  [C^(13)][l-C^]} .  A  rough  estimate,  similar  to  that  of  Eqs.  (23)- 
(28),  yields  something  like  exp(-4oL)  •  [4o!L  +  4 a£  )]  f°r  (30).  It 


30 


2 

shows,  using  (27),  that  the  terns  proportional  to  4a£  is  modified  dras- 
2 

tically  as  a ^  increases  above  unity.  Therefore  (28b)  is  not  valid, 
and  (20)  does  not  follow  in  this  case. 


The  following  approximation  is  the  basis  for  a  saturation-regime 

result:  Replace  all  Q  ,.^/k  by  the  function  p^z  )  to  be  defined  later, 

i  m 

The  index  i  pertains  to  the  B  factor  in  question.  Then  replace  C  (13) 

m 

by  an  average  <  0^(13)  >  .  All  M  factors  are  now  separable  and  we  can 
2 

reduce  <  I  >  to  an  exponential  form  in  the  same  way  that  (20)  was  de¬ 
rived  from  its  defining  series: 


<  l2  >  =  l2  exp{2(ke/4jt)2o/L  dz1/d2K]L®(K1)(l-C1 


) 


X  <  cos[K.Ap(z)]>+  0(o22)5/6) 


(31) 


It  is  reasonable  because  of  the  following  argument:,  £he  rough  estimate 


of  (30)  indicates  that  the  M-th  order  term  of  <  I4  >  is  of  order  (4a'L)M/M! 
2 

if  effects  of  0(a£  )  compared  to  CK L  are  ignored.  Terms  with  M  «  4Q!L 
(4o!L  is  a  very  large  number  compared  to  unity)  are  unimportant  and  con¬ 
sequent!),  Q_  is  a  sum  of  a  large  number  of  vectors  K. .  Furthermore, 

J?2^  ^^2^  ^  contains  a  sum  of  M-l  non -cancelling  vectors  (some  of 
these  occur  with  a  factor  2)  because  the  binary-correlation  property 
(17)  implies 

A(1V2)  ^0>-V4)  -  0  <32, 

for  all  m.  Consider  the  0^13)  factor  in  (30).  It  is  a  function  of 
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2  2 

of  ^(K^d  K2“— $(K^)d  K^j  with  the  other  M-l  filter  factors.  Because 
2-2C  (13)(1-C  )  «  2  in  the  integrand  it  follows  that  the  main  contri- 

m  m  “ 


m 


bution  of  C^(13)(l-C^)  is  found  by  netting  the  other  M-l  filter  factors 

equal  to  2  [see  the  discussion  immediately  following  (30)].  Consequently 

C^(13)  is  a  function  which  is  a  sum  of  all  vectors 

2  2 

of  the  remaining  integrand  ^(Kpd  K2 - ®(K  )d  Km>  There  is  complete 

analogy  to  a  two-dimensional  random-walk  problem  where  K.  is  a  random 

—  J 

phasor  and  ®(K.)  its  probability  density.  Therefore  Qri^-Q ^  tends 
to  a  Gaussian  random  variable  because  M  is  large,  and  with  zero  mean. 

We  replace  it  by  a  quantity  with  the  same  statistical  properties  and 
then  replace  C^(13)  by  <C^(13)>  .  This  Gaussian  quantity  will  be  named 
Ap(z)  at  present  and  it  will  be  identified  in  the  Appendix  with  ray  co¬ 
ordinates.  However,  at  present  it  need  not  be  discussed  further.  After 
having  made  the  substitution  forJ}2/k,  we  now  consider  the  M  =  2  factor 
analogous  to  (30),  and  repeat  the  very  same  argument  for  C2(13)  -  which 
involves  only  the  nK_)d  K„---tf’(K  )d  K  factors.  As  m  -+  M  it  appears 

incorrect  to  continue  this  procedure  because  the  number  M-l-m  of  K. 

—  j 

vectors  in  "^m+1^^  becomes  small.  Nevertheless,  the  error  in 

2 

so  doing  is  small  because  the  M-th  term  of  <  T.  >  is  considered 

together  with  many  subsequent  higher -order  terms  in  which  z  is  not  one 

m 

of  the  last  z  coordinates. 


N. 


A  straight-forward  extension  of  the  procedure  to  <1  >  yields 


<  IN>  =  IoNexp{2N(N-l)(keMjt)2  dz/c 
x<cos[K1-Ap(z1)]>+  N0(a22)5/6} 


jJd  K1#(K1)(1-C1) 


(33) 
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The  important  point  is  not  that^Q^/k  has  been  replaced  byjD(z) 

but  that  C  (13)  has  been  replaced  by  an  average  <  C  (13)  >  which  is 
m 

2 

not  a  function  ofJK^,--^.  If  we  may  ignore  0(c2  )  we  then  note  that 
we  have  shown  in  (33)  that  the  irradiance  is  also  log-normal  in  the 
saturation  regime. 


The  0(ct2  )  effects  become  negligible  as  the  Fresnel  number  L/kL^ 

goes  to  zero.  Consider  this  limit  for  the  sake  of  clarifying  (33). 

2  2 

As  -*•  “  (or  ->  oo)  the  irradiance  I  reduces  to  Iq  because  the  error 

2  1  /  6  2 

term  in  (33)  is  0(0^)“  '  .  The  only  way  to  allow  -*  00  consistent 
o 

wich  L/kL  0  and  wi;h  the  diffusion  condition  (12),  namely 
o 

1/3  2 

K  '  LC  «  1  is  to  allow  k  .  So  this  transition  is  essentially 
m  n  > 

to  a  high-frequency  limit  in  which  case  it  is  not  surprising  that 
I  ->  IQ  [17]. 

2 

However,  our  interest  is  in  the  case  that  k  =  const,  and  that  C 

n 

2 

^r.d/or  L  increase.  Let  us  take  L  ~  1  km  and  allow  to  increase  since 

this  is  in  accord  with  the  experimental  situation.  The  Fresnel  number 

24  2-42 

L/kLo  ~  10  at  optical  freqv  ncies.  Thus  a2  ~  10  and  the 

0(a2^)^^  terms  can  be  kept  very  small  in  (34)  compared  to  the  0(a£2)  *  ^ 

2 

terms  for  >10. 

Finally,  in  the  Appendix  we  work  out  an  estimate  of  (33).  It  yields, 


<  IN>  =  I  N  exp[2N(N-l)ac2-0(a12)"7/6] 
o  e  i 


The  last  result  indicates  that  the  log-amplitude  variance  decreases 
asymptotically  as  a  2,0(a.2)  ~  0(a  2)  We  discuss  this  further 


■  •  .  .  2 

in  section  V.  The  decrease  with  increasing  a  cannot  continue  indefinitely 

2 

because  ultimately  a 2  exceeds  unity  and  expression  (30)  is  no  longer 

2 

valid.  An  asymptotic  solution  for  the  case  that  -*  “  will  be  given  in 
the  next  section. 


IV.  THE  RAYLEIGH  REGIME:  k7^6LU^6C  2 

n 

Expressions  (24)  or,  equivalently,  (26)  have 
2 

differently  when  a2  »  1.  They  then  reduce  to 

taL[i1/2/27/3r(5/6)]-o21/3<!xp(-02) 


»  kL  2/L 
o 

to  be  handled  quite 


(35a) 

(35b) 


2 

Both  of  these  forms  vanish  as  -*  00  and  therefore  we  can  no  longer 
utilize  the  grouping  of  terms  in  (22).  Previously,  the  grouping  of  terms 
in  the  second  line  of  (22)  yielded  a  factor  2,  regardless  of  whether 
X)  =  0  or  AQ  t  0  in  any  of  the  C^ij).  Now  it  does  make  a  difference; 
for  AQ  0,  Cj(ij)  reduces  to  a  vanishingly  small  term  as  ->  in 
contrast  to  the  case  AQ  =  0  which  yields  C^ij)  =  1.  In  this  case,  we 
utilize  the  corollary  of  (18b)  which  differs  from  (18b)  by  an  exponential 
factor,  and  by  omission  of  i  =  j  factors.  Evidently,  the  main  contribu¬ 
tion  now  comes  from  a  subclass  of  diagrams:  those  for  which  all  C  (ill  =  1 

— —  mv  J  * 

i.e.  for  which  all  AQ  =  0.  Inspection  of  the  renormalized  beadless  dia¬ 
grams  in  conjunction  with  the  property  (32)  shows  that  this  subclass  of 
diagrams  has  the  following  properties: 
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(i)  There  are  no  <  13 >  nor  <  24  >  rungs  (the  notation  <  13  >  in¬ 
dicates  correlation  between  the  first  and  the  third  horizontal 
axes).  In  general,  there  are'  no  <  BB  >  nor  <  B*B*  >  rungs. 

(ii)  Any  B  axis  is  connected  by  rungs  to  only  one  B*  axis,  and 
vice  versa.  That  is  to  say,  all  diagrams  can  be  factored 
into  products  of  <  BB*  >  diagrams . 


Instead  of  (22)  and  similar  factors  for  m  =  2,  —  ,M  as  the  contri- 

2 

bution  of  a  diagram  with  M  rungs  to  <  I  >  ,  we  have  in  general  for 
N 

<1  >  according  to  the  above  procedure, 


M  j2  2  rra-1 
exp(-2NQ!L)  •  IT  - _  „ 

m=1  16jf2  ° 


r\  f< 


dz  Jd  K  ®(K  )C  (ij) 
m  m  mm 


M, 


=  exp(-2NC&i)  *(2o!L)  /Ml 


(36) 


Let  us  try  to  add  up  (36)  for  all  M-rung  diagrams  where  n  pairs  of  B 
and  B*  axes  have  no  rungs.  We  may  not  have  any  mixing,  i.e.,  we  cannot 
have  two  rungs  connecting  one  B  to  two  B*  axes,  or  vice  versa,  otherwise 
at  least  one  Cm(ij)  ^  1.  For  one  choice  of  n  pairs  of  conjugate  B  and 
B*  axes  there  are  (N-n)!  ways  of  pairing  the  remaining  N-n  pairs  of  B 
and  B*  axes  two  by  two.  Let  us  assume  for  one  choice  of  pairing  that 
there  are  m^n^, - "'N-n  run8s  on  N-n  pairs  we  have  chosen.  Ob¬ 
viously  m^  +  —  +  =  M.  Then  we  obtain 
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(37) 


(N-n)'. 


e -2 Noil  f  ___  £  M! 

V*  mN-n=1  mr  —  mN-n' 
=  (N-n)'.  e2NaL(l-e'20iL)N_n  , 


-MQL 

e 


for  all  M-rung  diagrams  that  leave  a  single  choice  of  n  pairs  of  B,B* 
axes  unconnected  to  each  other  by  rungs  >  We  need  one  more  numerical 
factor  in  (37)  before  summing  over  n  from  n  =  0  to  n  =  N,  namely  the 
number  of  ways  that  we  can  choose  n  pairs  B  and  B*  from  N  pairs  B  and 
B*.  This  number  is  easily  seen  to  be 


2  2 
NZ(N-1T 


2  2 

lV 


---  (N-n-fl)2 
2 

--  n 


(N!)2 

(n!)2[(N-n)!]2 


(38) 


Thus,  we  obtain  <  I^>/<I>N  by  summing  the  product  of  (37)  with  (38) 
over  n  from  n  =  0  to  n  =  N: 


<IN>  =  <  I>N 


f.  -,.-(N'.)2 

(n'.)2(N-n)! 


-2nQ!L  -20L.N-n 
e  (1-e  ) 


(39) 


The  meaning  of  (39)  for  small  and  intermediate  OIL  has  been  worked  out 
in  [2]  and  [3],  but  in  the  present  case  (XL  »  1  and  we  retain  only  the 
n  =  0  term  in  (39)  to  obtain, 


<IN>  =  N'.  IoN[l  +  0(a2l/3eP2)] 


(40) 


This  result  is  the  Rayleigh  limit;  a  special  case  of  a  result  found 
ear^er  ^or  region  IV  of  Fig.  1  where  (XL  had  not  greatly  exceeded  unity. 
It  states  that  as  a2  00 »  the  irradiance  approaches  an  exponential  dis¬ 
tribution  (i.e.,  the  amplitude  becomes  Rayleigh-distributed). 


36 


V.  COMMENTS  ON  THE  IRRADIANCE -VARIANCE  SATURATION 

The  results  (3)  -  or  more  fully  (21),  (28),  (29),  (34),  and  (40)  - 

allow  us  to  reconsider  the  ir radiance -variance  saturation.  It  seems 

somewhat  easier  to  consider  the  log-amplitude  X  =  %  In  I  and  to  plot 
2  2 

<X  >  as  a  function  of  .  The  results  indicate: 

(i)  <X^>  =  a ^  for  values  of  a ^  much  less  than  (fc  ^L/k)  t 
a  number  of  order  unity. 

(ii)  <  X^>  «  o  ^'(o^)  ^or  intermediate  values 

(K^L/k)  «  c ^  «  kL^/L.  The  log-amplitude  variance 

2  ••1/6 

changes  from  linear  growth  to  a  decrease  as  (a ^  )  '  as 

2 

a£  grows  beyond  a  number  of  order  unity.  This  trend  appears 
to  agree  well  with  experimental  observation. 

.  .  2 

(iii)  <  X  >  should  level  off  asymptotically  at  the  level 
2 

x  /24  «  0.41  when  the  Rayleigh  limit  is  reached,  i.e.,  when 

2  2  2 
»  kLQ  /L.  It  depends  on  the  magnitude  of  kLQ  /L 

2  2 

whether  <  X  >  has  dropped  below  jt  /24  before  leveling  off 
or  not . 

Experimental  data  now  exists  with  which  to  check  these  aspects  of 
the  theory.  Mevers,  et  al.  [18]  have  found  a  decrease  with  (a^) 
in  region  (ii),  in  good  agreement  with  our  result.  The  newer  data  of 
Kleen  and  Ochs  [19]  as  well  as  the  older  data  of  Gracheva  [20]  exhibit 
a  trend  of  this  type  (although  no  analytical  expression  for  their  data 
is  given  and  although  their  data  exhibit  appreciable  spread).  Furthermore, 
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L  2 

<X  >  appears  to  drop  below  the  ultimate  level  of  it  /24  (presumably 
2  ^ 

because  kL^  /L  ~  10  i.s  so  large  a  number) . 


Figure  4  gives  a  sketch  of  a  saturation  curve  satisfying  the  above 

2 

predictions.  The  dip  below  the  Rayleigh  limit  it  ./24  has  been  shown  in 


Figure  4.  Log-amplitude  variance  <y?>  vs  a ^  %  \J ^  Cn^: 

A  sample  curve  fitting  the  (dashedf  asymptotes  given 
by  theory. 


exaggeration  [it  may  or  may  not  occur,  depending  upon  the  magnitude  of 

2  "1/6 

the  coefficient  of  (ce  )  ].  The  major  weakness  of  the  theory  is  its 
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inability  to  predict  the  coefficient 

2 

ration  limit.  The  maximum  of  <  X  > 
in  the  present  form  of  this  work. 


of  (a/)’176  in  <  X2> 
is  therefore  also  not 


in  the  satu- 
predictable 
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APPENDIX:  SATURATION -REGIME  LOG -AMPLITUDE  VARIANCE 


The  result  of  section  III  for  the  saturation  regime  was, 

<  IN>=  IoNexp[2N(N-l)<SX2>  ], 

2<SX2>=  (ke/4it)2  J*  dzJd2K®(K){l-cos[K2(L-z)/k]}  (Al) 

x  <  cos[K*Ap(z)]>  . 

0m  mm 

Here,  it  will  be  argued  that  Ap(z)  is  the  difference  of  two  ray-coordin- 

2 

ate  realizations,  and  a  first-order  estimate  of  <  5X  >  will  be  given. 

In  order  to  determine  the  nature  of  Ap(z)  we  will  rederive  (Al)  in  another 
fashion  by  regarding  the  equation  for  i|r  =  InB  obtained  from  the  parabolic 
wave  equation, 

2ik  ty/dz  +  +  k2Se  +  (VTt)2  =  0  .  (A2) 

The  Rytov  approximation  (2)  results  when  the  last  term  of  (A2)  is  dropped. 
We  will  follow  a  procedure  analogous  to  that  described  in  [4]  to  circum¬ 
vent  the  difficulty  of  the  last  term.  Let  \|r  =  X  +  icp.  We  obtain  the 
ray  equation  kdp(z)/dz  =  from  (4)  by  modifying  the  stationary-phase 
method  to  the  more  precise  steepest-descent  approximation  [4].  Consider 
the  transformation  from  the  cartesian  coordinate  system  (x#y,z)  to  the 
locally  orthogonal  curvilinear  system  (p(s),s).  It  is  easily  seen  that 

mmrn 

=  0  in  the  new  coordinate  system;  'v5j,i|r  and  are  invariant  with  re¬ 
spect  to  the  transformation  and  they  may  be  expressed  in  any  coordinate 
system.  However  cfy/dz  is  not  invariant.  Nevertheless,  the  curvilinear 
effects  can  be  ignored  in  cfy/dz  because  the  difference  between  d\jr/dz  and 
cty/ds  is  LC^  )  as  given  by  (11).  Thus  to  good  approximation 
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(A3) 


2ik  cfy/ds  +  +  k2fie  =  0 

along  rays  determined  by  kdp(z)/dz  =  in  terms  of  the  original  coordin- 
ate  system.  By  utilizing  fie  (s)  =  fie(p(z),z)  and  by  replacing  dz  «  ds 
and  (L-s)  «  (L-z)  we  write  the  solution  of  (A3)  as, 

\| t  =  (ik/Sn2)^  dz/d£(K,z)exp[-iK*[p(z)+(L-z)K/2k]}  .  (A4) 

Eq.  (A4)  constitutes  an  important  correction  to  (2),  even  in  the  limit 

p  -*■  0,  because  <\| t>  t  0.  Tatarski,  in  a  slightly  different  formulation 

[21]  has  computed  <X>  in  the  geometrical-optics  regime  where 
2  2 

sin[K  (L-z)/2k]  may  be  replaced  by  k  (L-z)/2k.  It  is  this  difference 

2 

with  (2)  that  yields  the  first  equation  in  (Al),  where  <&X  >  repre¬ 
sents  the  variance  of  X  =  Rety.  The  derivation  of  (34)  ensures  the  con- 

2 

servation  law  <  X  >  +  <8X  >  =  0  so  that  it  suffices  to  compute  one  of 

2 

these  two  quantities.  If  <  SX  >  is  formed  from  (A4)  by  subtracting 
2  2 

<X>  from  <X  >  ,  we  obtain  (Al)  because  the  correlation  between 
d£(K1,z1)  and  p(z_)  can  be  ignored  (note  that  the  multivariate  Gaussian, 
or  Markov,  property  holds  for  these  quantities,  e.g.,  as  in  paragraph 
65  of  [16]).  The  modif ied-Rytov  approximation  (21)  follows  by  setting 
Ap(z)  =  0 . 

The  above  procedure  clearly  is  an  improvement  upon  the  customary 
Rytov  approximation.  However  (A4)  is  valid  beyond  the  region  in  which 
Cj  «  1.  It  breaks  down  when  C2  ,  given  by  (29),  approaches  unity  be¬ 
cause  the  WKB  expression  for  phase  is  significantly  different  from  the 
straightened -ray  expression.  This  can  be  seen  from  (A4)  by  ignoring 
the  filter  factor  and  repeating  the  estimate  of  (23a)  with  £QL/k  replaced 
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by  Ap  (with  variance  ~  ^1/3L3Cn2).  Thus,  accepting  the  validity  of 
(A4)  for  (7*2  «  1,  we  find  (Al)  from  (A4)  as  described  above.  By 
comparing  the  result  of  this  method  to  (33)  we  have,  been  able  to  identify 
the  hitherto  unspecified  Ap(z)  with  the  difference  of  two  ray-coordinate 
realizations . 

Now  an  estimate  of  <SX2>  can  be  given.  In  (Al)  we  set 
<cos[K*Ap(z)]  >  =  exp[-K2  <Ap2(z)>  /4]  by  utilizing  the  Gaussian  property 
of  A£.  Furthermore,  we  insert  e2«(K)  «  32jtrCn2K-11^3  from  (25)  into  (Al), 
transform  variable  K  into  x  =  K2(L-z)/k,  to  obtain 

<SX  >  =  rk7/6C^  tffLdzz5/6o/°dxx-11/6(l-cosx)exp[-xk  <Ap2(z)>/4(L-z)]  , 

(A5) 

where  coefficient  y  83  0.033*2.  The  major  problem  is  to  compute  <  Ap2(z) >  . 

The  ray  coordinates  ^p(z)  converge  to^(L)  =  0  and  diverge  randomly  as 

L-z  grows.  Considering  the  fact  that  the  integrand  of  (Al)  is  negligible 

for  L-z  -  0,  we  shall  approximate  <Ap2(z)>  «  4<p2(z)>  as  if  both 

ray  realizations  are  independent.  This  approximation  should  be  better 

2 

for  large  L,  i.e.  for  large  .  We  utilize  the  well-known  geometrical- 
optics  expression  for  <p2(z)>  given  by  Tatarski  [21]  when  p(L)  =  0 
and  ray  derivative  p’(0)  =  0, 

<p2(z)>  =  (4r/3)Cn2/Cm1/3[(L-z)3  +  3z(L-z)2]  (A6) 

to  obtain, 

<5X2>  =  rCn2k7/6L11/6o/dyy5«/te-ll/6(1.co3j[)exp[^i2f(y)] 
f(y)  =  (4y/3)(l-y)(l+2y)  ,  (A7) 
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221/32  2 

and  0.  =  kL  K  C  in  agreement  with  (29).  It  can  be  seen  that  <5X  > 

i  m  n 

reduces  to  a  2  =  0.31k^^L^^^C  2  in  the  limit  a  2  -*  0.  In  order  to 
e  n  1 

2 

obtain  an  expression  for  large  ,  we  return  to  (Al)  and  partially  in¬ 
tegrate  (....)  with  respect  to  z.  Upon  utilizing  the  approximations 
incorporated  into  (A6)  and  (A7)  we  obtain 


<  SX2>  =  rCn2k7/6L11/6(  f  dxx-11/6(l-x'1sinx)exp(-  |  ya^x) 

o 

(A8) 

+  8i raj2  ^  dyy(l-y)  j  dxx  5^6[l-  Jexpt-xa^fCy)]) 


2 

For  large  ,  it  is  permissible  to  replace  the  filter  factors  in  each 

2 

term  of  (A8)  by  the  x  term  of  their  power-series  expansion.  It  can  be 

seen  from  the  ensv.ing  expression  that  the  leading  terms  of  (A8)  for  large 
2 

Oj  yield  the  dependence. 


2 

SX  >  = 


•0<OlV7'6 


which  we  have  utilized  in  (34).  The  numerical  coefficient  can  be  com¬ 
puted  from  (A8),  but  it  may  only  have  an  approximate  meaning  in  view 

2  2 

of  the  approximation  <  Ap  (z)>  «  4<p  (z)>  . 


Note:  A  revised  version  of  pp.  9-43  is  currently  being  prepared 
for  publication. 
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V.  I.  Klyatskin  in  a  recent  paper  [Sov.  Phys.-JETP  33,  703  (1971)] 

N  N 

has  come  to  the  different  conclusion  that  <  I  >  -*  N<I>  as 
2 

-*  and  that  no  probability  distribution  exists  in  this  limit. 

This  conclusion  appears  not  valid  for  optical  wave  propagation  in 

2 

turbulent  air  for  the  following  reasons:  The  limit  ag  -*  00  violates 

1/3  2 

either  the  no -diffusion  approximation  Km  LC^  «  1  implicit  in 

the  parabolic  wave  equation,  or  the  smallness  of  the  macroscale 

2 

Fresnel  number  L/kLo  ,  or  it  corresponds  to  the  trivial  case 

k  -*•«.  Klyatskin1  s  treatment  ignores  the  effect  of  macroscale 

L  ,  hence  it  deals  with  an  inf in ite -macro scale  medium.  In  that 
o 

case  it  is  not  surprising  that  the  Rayleigh  limit  (see  Section 

IV)  is  not  found.  For  the  then  non-trivial  case  that  k  is  finite 

Klyatskin' s  result  can  only  hold  for  an  infinite-macroscale  medium 
2  2  7 /6 

in  which  Cn  L  is  finite  but  L  '  -*•  »;  a  clear  impossibility 

2 

because  must  be  independent  of  L.  If  on  the  other  hand,  we 

allow  k  -»  a>,  it  is  physically  obvious  that  the  free -space  irradi- 

ance  I  must  result . 
o 
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